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Abstract 

Hadronic axions with the decay constant fa — 10^ GeV may fulfill all 
astrophysical and laboratory constraints discussed so far. In this paper, we 
reexamine the possibility of the hadronic axion window while taking into 
account the uncertainties of some parameters describing low energy axion 
dynamics. It is found that /„ in the range from 3 x 10^ GeV to 3 x 10^ can 
not be excluded by existing arguments. We then examine the implication 
of this hadronic axion window for the big-bang nucleosynthesis (NS) by 
evaluating the energy density of thermal axions at the nucleosynthesis 
epoch. Our analysis yields {pa/ Pv)ns = 0.4 ~ 0.5 which exceeds slightly 
the current best bound {pa/Pu)NS < 0.3. 



1. Introduction 

One of the most attractive solutions to the strong CP problem is to introduce 
a spontaneously broken Peccei-Quinn symmetry which gives 6 = dynamically 
This solution predicts the existence of a pseudo-Goldstone boson, the axion, whose 
decay constant fa is tightly constrained by astrophysical and cosmological consider- 
ations Pi . One frequently quoted window for the axion decay constant satisfying all 
phenomenological constraints is = lO^*' ~ 10^^ GeV. Besides this, for hadronic 
axions which have vanishing tree level coupling to the electron, there can be another 
window (the hadronic axion window) around fa — 10® GeV. The existence of the 
hadronic axion window relies upon (i) the axion-electron coupling is radiatively in- 
duced and thus highly suppressed ||^, (ii) axions on the window are trapped inside the 
supernova core 0], and (iii) a significant cancellation can occur in the axion-photon 
coupling between the short distance contribution (from the electromagnetic anomaly 
at fa) and the long distance contribution (from the axion-pion mixing at the QCD 
scale) Based upon these points, fa — 10® GeV can be consistent with a variety 
of the strong astrophysical constraints including those from the supernova SN1987A 

i- 

The hadronic axion window can be relevant for some models of dynamical sym- 
metry breaking. Then it would become important to examine the phenomenological 
viability of the window more carefully. To this end, recently the effect of relic hadronic 
axions upon the diffuse extragalactic background radiation has been studied exten- 
sively 1^. By examining the photon flux from decaying axions 0, it has been noted 
that a severe cancellation should occur in axion-photon coupling for fa — 10® GeV to 
be consistent with the observed data |^. In fact, there is a possible loophole for this 
conclusion. Even when fa and the axion-photon coupling are fixed, there still exists a 
large uncertainty of the photon flux associated with the uncertainty of the light quark 
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mass ratio z = iriu/md. For a given value of /„, the axion mass is determined by z (see 
eq. (5)). In ref. 0, z = 0.56 was taken to arrive at the above conclusion. However it 
has been argued by several authors that, due to instanton effects, the true value 
of z can be significantly smaller than the conventionally used value 0.56. A smaller 
z would give a smaller axion mass for which both the position and the intensity of 
the photon flux become lower. Clearly then the required cancellation of the axion- 
photon coupling becomes weaken. At any rate, in generic hadronic axion models, one 
can achieve the required cancellation by appropriately arranging the electromagnetic 
anomaly at fa- 
in fact, the axion-photon coupling is the least known parameter among those 
describing the low energy dynamics of hadronic axions. It is thus desirable to consider 
an effect which is independent of the axion-photon coupling. In this paper, we wish 
to examine one such effect, that on the big-bang nucleosynthesis^. The organization 
of this paper is as follows. In the next section, we briefly review low energy axion 
couplings [|, |] which will be used in later discussion. In sec 3 we collect arguments 
defining the hadronic axion window. Taking into account the uncertainties of the 
involved parameters, we argue that existing arguments allow the window: 3 x 10^ < 
/a < 3 X 10^ (in GeV unit). In sec 4 we evaluate the rate of thermal axion production 
after the quark-hadron phase transition which was roughly estimated in ref. ||10| . This 



will determine the energy density of thermal axions at the nucleosynthesis epoch. For 
the window defined in sec 3, we find {pa/Pv)NS = 0.4 ~ 0.5 which exceeds shghtly 
the current best bound {pa/ Pu) < 0.3 ||lT|. Sec 5 is given for the conclusion. 

2. Low energy couplings of hadronic axion 

In this section, we review the low energy couplings of hadronic axions P, ^ . We 
^ The constraints imposed by the nucleosynthesis on the interactions of majorons have been 
discussed by Bertohni and Steigman B 
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will mainly follow ref. [H. To be definite, let us consider the simplest hadronic axion 



model 1 12 1 which was first considered by Kim. The model contains an isosinglet Dirac 
quark Q with the electromagnetic charge Yq and also a gauge singlet complex scalar 
field a. These exotic fields carry nonzero PQ charges: 1, —1, and 2 for Ql, and 
a respectively, while all other fields have vanishing PQ charge. The PQ symmetry is 
spontaneously broken by (a) = fa/ \/2 and then the axion a can be defined by 

= ^{fa + p)exp{ia/fa). (1) 

Integrating out heavy degrees of freedom, viz Q and p whose masses are of the order 
of fa, we obtain the axion effective lagrangian at the renormalization point ^ = fa'- 



1 2 2 

2^ ^ ' /a 327r2 327r2 ' ^ ' 



Here G and Y denote the field strengths of SU{?>)c and f/(l)y respectively, while G 
and Y are their duals, and n = QYq denotes the electromagnetic anomaly arising from 
the triangle loop of Q. Note that there is no tree level axion-electron coupling, which 
is the characteristic property of hadronic axions. 

Given the axion effective lagrangian at /i = /a, we can derive a variety of axion 
properties at low energies by scaling the effective lagrangian down to the QCD scale 
A ~ 1 GeV 0. Below A, the QCD anomaly and instanton effects become strong and 
also the chiral symmetry is spontaneously broken. It is then convenient to eliminate 
the aGG coupling through an axion- dependent chiral rotation of the light quark fields 
q = (m, d, s): 

qL ^ exp{iaQA/fa)qL, qn ^ exp{-iaQA/ fa)qR, (3) 

where 

Qa = M-V2tr(M-i), (4) 
and M = diag(m„, m^, m^) denotes the light quark mass matrix. With this chiral 



rotation, the axion field a becomes a physical mass eigenstate with 



""""1 + . fa ' 

where = 93 MeV denotes the pion decay constant and z = rriu/md. Also the axion 
couplings with the electron and the photon renormalized at /i < A are given by 

Ca^^TT^F'^^F,, + C^^^eri.e, (6) 

Ja 327r^ fa 

where0 

2(4 + z) 



Cn/y — K, 



3(1 + z) 

a. = ^[«:ln(^) - l^ln(-)]. (7) 

There are two kind of axion couplings with hadrons; derivative and nonderivative 
couplings In fact, axion phenomenology associated with the hadronic axion win- 
dow is governed mainly by derivative couplings. At /i ~ A, after the chiral rotation 
(3), axion derivative couplings with the light quarks are given by 

-(1 + 7)g7^75g + 2^ag7^75A'^g], (8) 

where Xa = tr((5AA") and 7(/i) is due to the renormalization of the singlet current 
occurring between fa and /i. Here tr(A'*A^) = 25ah- This then gives 

^(^^o" + ^a;.J,-), (9) 

where Jq and are the hadronic currents made of the pseudoscalar meson octet 
S = exp(i7r'^A"//^) and the baryon octet B = B'^X^'/V^: 

= ^/2tr(A"(SD^St-StD^S))+tr(57^[A?,(0,5]) 
+FtT{BYl5[XsiO,B]) + DtT{BYldm),B}), 
= Sti{Bri,B). (10) 



Note that our result for Cae is smaller than the Srednicki's in ref. H by the factor 1/2it. 



Here X%{0 = (^A'^C^ - e^A«0/2, X^O = (^A^^^ + e^A«0/2 for e = S, and F and 
denote the baryon matrix elements of the SU (3) octet axial vector currents, while 5" 
is defined as 

(1 + 7) < B\q-i^'-ir,q\B >= Sn", (11) 

where denotes the covariant spin vector of the baryon state B. Note that S is 
independent of the QCD renormalization point 

Among axion couplings with hadrons, those of particular interest for us arc the 
axion-nucleon couplings CaN responsible for the processes NN NNa and Ntt — > 
Na and also the axion-pion couplings CaTv for tttt — > tto. These couplings can be read 
from eqs. (9) and (10) giving 

J a J TV 

jCa. = C„,|^(7r°7r+9^7r- + 7rV-9^7r+-27r+7r-9^7r°) (12) 

JaJiT 



where 



^an — T^TT- + ^TT" T-D + -S, 

2(1 + 2;) 6(1 + 2;) 6 



Cap = — ^— F + 4 — '^D + -S, 
2{l + z) 6(1 + ^) 6 ' 



2V2{l + z) 

^"'^ ^ 3(TT^ ^^^^ 

Several parameters appear in eqs. (7) and (13) to determine the low energy axion 
couphngs for a given value of /„. Let us discuss the values of those parameters. First 
of all, the electromagnetic anomaly k is due to the physics at the axion scale /„ and 
thus is totally unknown. In principle, one can arrange the particle content at /„ to 
have K taking an arbitrary value of order unity. The axial vector coupling constants 
F and D of baryons can be determined by the nucleon and hyperon beta decays 
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as F = 0.47 and D = 0.81 [|T3|. Then recent EMC measurement^ of the polarized 
proton structure function gives S* ~ 0.13 ± 0.2 |1^. About z, one usually uses the 
result z = 0.56 from first order chiral perturbation theory |T^. In second order chiral 



perturbation theory, z receives a correction 6z = 0{ms/ATTfT^) which is essentially due 
to instanton effects 0. A larger value of this instanton- induced correction means a 
smaller z. It has been argued that 6z can be large enough to imply z = 0. With 
this point taken into account, we will allow z take a value significantly smaller than 
0.56. 

3. The window 

In this section, we will collect arguments which define the hadronic axion window. 
Taking into account the uncertainties of k and z, we will argue that existing arguments 
do not exclude the window 

3 X 10^ GeV < fa <Sx 10^ GeV. (14) 

A key property of hadronic axion is that its coupling to the electron is radiatively 
induced and thus is highly suppressed. As is well known, helium ignition of red 



giants provides a strong constraint on the axion-electron coupling |jT^, C*ae^e//a < 
1.5 X 10^^'^ in our notation. This gives^ fa^KX 10^ GeV which is satisfied by the 
window (14) for n of order unity. 

The list of other arguments defining the hadronic axion window is as follows: 



^ In fact, what the EMC measurement gives is the value of 5/(1 + 7) at /i ~ 3 GeV where 7 is 

negligibly small. Note that at leading order the evolution of 7 is given by fidj/dii ~ — 3(ac/7r)'^7 

[ p^ with j{fi = fa) ~ 0. It has been pointed out Q that elastic neutral current experiment gives 

a similar result S ~ 0.09 ± 0.24. 

* Since our result for the radiatively induced axion-electron coupling Cae is smaller than the 

Srednicki's, the lower bound on fa is relaxed by the factor 1/2tt. Note that the hadronic axion 

window becomes widen as a result of this point. 
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(A) the axion's effect on the neutrino burst from SN1987A [T^, (B) the effect of 
axions emitted from SN1987A on the Kamiokande II detector [jT9|, (C) the effects 
of axion emission on the evolution of hehum burning low-mass stars [^], (D) the 
effect of decaying relic axions on the diffuse extragalactic background radiation 
In regard to the argument (A), axions on the window (14) would be trapped inside 
the supernova core @]. A detailed numerical analysis for this "trapping regime" has 
been performed by Burrows et al. JlSl who have concluded that if fa < 8CaN x 10^ 
GeV (here CaN — (C^n + ^%}}^^'^ I ^) axions would be so strongly trapped that 
axion emission would not have a significant effect on the neutrino burst[| Also by 
considering (B), it has been pointed out that < ICaN x 10^ GeV would have 



produced an unacceptably large signal at the Kamiokande detector [T^. These then 
give the allowed window^ 

2CaN X 10^ GeV < /a < 8Ca7v X 10^ GeV. (15) 

Using F = 0.47, D = 0.81, and S = 0.13 ± 0.2, and allowing z vary between zero 
and 0.56, we find CaN = 0.15 ~ 0.4 for which the above window goes into that of eq. 
(14). 

Can the window (14) be compatible with the constraints from (C) and (D)? The 
arguments (C) and (D) provide limits to the axion-photon coupling. From (C), one 
has derived |^ < (/^/lO^GeV). It has been pointed out that, for z = rriu/md = 
0.56 and fa — 10^ GeV, the argument (D) can lead to a more stringent limit ||^, 
say Ca-y < 0.02. However the validity of the latter limit strongly relies upon the 
used value of z. Since rria oc ^/z, if \z\ <^ 0.56 which is an open possibility |^, 
both the position and the intensity of the photon-line from decaying axions would 



^ Although less complete, axion trapping was discussed also by Ishizuka and Yoshimura pi| and 
Engel et al. 

^Here we ignore possible uncertainties of the astrophysical arguments leading to this range of fa- 
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be significantly lowered. Clearly then the limit from (D) becomes invalid or at least 
weaken. 

For hadronic axions, can be divided into two pieces, the short distance con- 
tribution K denoting the electromagnetic anomaly at the axion scale fa and the long 
distance part —2(4 + z)/'i{l + z) arising from the axion-pion mixing. Then an inter- 
esting point is there can be a severe cancellation between the short and long distance 
contributions. In order to see this, let us consider two examples: (e.l) z = 0.56, 
Yq = 4/7; (e.2) z = 0.05, Yq = 2/3. Here Yq denotes the electromagnetic charge 
of the heavy Dirac quark Q and thus k = 6Yq. We then have (e.l) Ca-y = O(10~^); 
(e.2) Ca-y ^0.1. These examples simply show that a significant cancellation can occur 
rather naturally, but the degree of cancellation is highly model-dependent. At any 
rate, it is true that, for on the window (14), one can always arrange the model to 
have K for which Ca-y becomes small enough to satisfy the constraints from (C) and 
(D). We thus conclude that, due to our ignorance of k ( and also of z), the hadronic 
axion window (14) is allowed also by (C) and (D). 

4. Thermal hadronic axions and the nucleosynthesis 

In this section, we consider the implication of the hadronic axion window (14) for 
the big-bang nucleosynthesis (NS). As is well known, in the standard model of NS, 
the energy density of exotic particles at the NS epoch is severely constrained. The 
primordial yields of ^He is sensitive to the universal expansion rate H oc p^/^ at T ~ 1 
MeV and, thereby, leads to a bound on any exotic contribution to the energy density 
p. The recent work of Walker et al. |jll[ provides the most stringent bound on the 
axion energy density normalized to the energy density of the electron-neutrino: 

6N^ = ( ^ ) < 0.3. (16) 

Since both axions and neutrinos are relativistic at T ~ 1 MeV, one simply has 
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SNy = {A/7){Ta/TyYNg. Suppose axions were initially in equilibrium with the thermal 
bath of normal particles (7, e^, i/, ...) at high temperature, but later become decouple 
at temperature To > 1 MeV. Entropy conservation then leads to the relation 



where (7* counts interacting degrees of freedom whose entropy density is given by s = 
27r^(7*Tj^/45. Here the axion decoupling temperature Td can be found by comparing 
the expansion rate H with the axion production rate 



where denotes the axion number density and the sum extends to all production 
processes involving as initial states the particles i and j which are still in equilibrium 
at Td > 1 MeV. As was pointed out by Bertolini and Steigman for the NS 
constraint (16) to be satisfied, axions (or any light scalar with m ^ 1 MeV) must 
decouple before T = 100 MeV. 

Clearly ^A*"^ < 4/7. If axions decouple before the quark-hadron phase transition, 
we would have ^A*"^ < 0.06 because of the relative "heating" of neutrinos during 
the quark-hadron phase transition. We are thus led to consider the possibility for 
axions in thermal equilibrium after the quark-hadron phase transition. Then there 
are two types of processes dominantly producing thermal axions: (a) tt^tt^ ^ ott^, 
TT+TT" — i> avr", and (b) ti^N aN {N = n or p), vr+n ap, 7r~p — >■ an and their 
CP conjugates. In fact, the rate of the process (b) was roughly estimated in ref. @]. 
Here we will carefully evaluate the rates of both (a) and (b). 

Let T(^a) and r(ft) denote the axion production rate of (a) and (b) respectively. 
Then a straightforward computation yields 




(17) 



T = n 



-1 



(18) 



'a 




3 1 
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= ^y^[^7i(5CL + 5Ci + 2anC^ap) + 6CL^]/(,), (19) 

where 

/(,) = n-'T' / dx,dx2^f{yr)f{y2) / t/o;^ ^^^^-^ ^, 

J yiy2 J-i s^T* 

h) = n-^n^T'^ j dxixivlfiyi). (20) 

Here f{y) = l/(e^ — 1) denotes the pion distribution function, Ua and are the 
axion and nucleon number densities at equihbrium, Xi = \pi\/T, yi = Ei/T {i = 1,2), 
and s = 2(m^ + T^(?/i?/2 — xiX2Uj)). To derive T(^a}, we have used the pion- nucleon 
couphng QAN'^^'^^d^TxN where qa = F + D and vr = 7rVY2 (r* = the Pauh matrices). 

By performing numerical analysis for /(a) and /(&), one can see that T{^a) dominates 
over Tqj^ when T < 150 MeV. Here to be definite, we have used z = 0.56, F = 0.47, 
D = 0.81, and S = 0.13, but the result is quite insensitive to the allowed variations 
of these parameters. Using r/i7~latT = To, we findQ To = 30 ~ 50 MeV for the 
range of fa from 3 x 10^ GeV to 3 x 10^ GeV. With the relation (17), this gives 

= 0.4 ~ 0.5, (21) 

which is in slight contradiction to the NS bound (16). 

5. Conclusion 

Hadronic axions with the decay constant around 10® GeV can escape from all 
astrophysical and laboratory constraints discussed so far. In this paper, we have col- 
lected arguments which define the hadronic axion window and found that existing 
^ Using the rough estimate of Fj-;,) given in ref. Td ~ 50 MeV was found in ref. Our 

careful evaluation yields Ff^) which is significantly smaller than that of ref. This is mainly due 
to the suppression caused by the axion-nucleon coupling constants. As a result, if one uses our F(f,) 
while ignoring T(^a), one would get To = 70 - 100 MeV for /a = 3 x lO'^ - 3 x 10^ GeV. However 
as we have mentioned, around Td, the pion process (a) largely dominates over the nucleon process 
(b), leading to a lower To. 
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arguments allow the window 3 x 10^ < /a < 3 x 10^ (in GeV unit). The phenomeno- 
logical viability of this window relies strongly upon the mo del- dependent cancellation 
that occurs in the axion-photon coupling. As we have noted, in generic hadronic 
axion models, it is not so unnatural to achieve the required cancellation. Since the 
axion-photon coupling is quite model-dependent and is totally unknown, it is desir- 
able to consider an effect which is independent of the axion-photon coupling. We 
thus have considered the effect on the nucleosynthesis which is determined by the 
axion-hadron couplings. By evaluating the axion production rate after the quark- 
hadron phase transition, we could determine the energy density of thermal axions 
at the nucleosynthesis epoch. It is quite certain that To < 100 MeV and thus the 



resulting 5Ny exceeds the lower bound 0.3 |Tl|]. However since the deviation is not so 
significant, 5N,j = 0.4 ~ 0.5, it would still be difficult to conclude that the hadronic 
axion window is excluded by the primordial abundance of ^He. 
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